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Abstract 



^ ■ We numerically reexamine the scaling behavior of period doublings in four- 

dimensional volume-preserving maps in order to resolve a discrepancy between 
numerical results on scaling of the coupling parameter and the approximate 
On . renormalization results reported by Mao and Greene [Phys. Rev. A 35, 3911 

J>^. (1987)]. In order to see the fine structure of period doublings, we extend 

^ . 
Q . the simple one-term scaling law to a two-term scaling law. Thus we find 

,^ . a new scaling factor associated with coupling and confirm the approximate 

renormalization results. 
X ■ PACS numbers: 05.45. +b, 03.20.+i, 05.70.Jk 
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Universal scaling behavior of period doubling has been found in area-preserving maps 
||T|-^. As a nonlinearity parameter is varied, an initially stable periodic orbit may lose its 
stability and give rise to the birth of a stable period-doubled orbit. An infinite sequence of 
such bifurcations accumulates at a finite parameter value and exhibits a universal limiting 
behavior. However these limiting scaling behaviors are different from those for the one- 
dimensional dissipative case P]. 

An interesting question is whether the scaling results of area-preserving maps carry over 
higher-dimensional volume-preserving maps. Thus period doubling in four-dimensional (4D) 
volume-preserving maps has been much studied in recent years [|7|,|9Hl3|. It has been found 



in Refs. ||TTHl3[ that the critical scaling behaviors of period doublings for two symmetrically 



coupled area-preserving maps are much richer than those for the uncoupled area-preserving 
case. There exist an infinite number of critical points in the space of the nonlinearity and 



coupling parameters. It has been numerically found in [|ri|)0 that the critical behaviors 
at those critical points are characterized by two scaling factors, 6i and 62- The value of 61 
associated with scaling of the nonlinearity parameter is always the same as that of the scaling 
factor 6 (= 8.721 . . .) for the area-preserving maps. However the values of 62 associated with 
scaling of the coupling parameter vary depending on the type of bifurcation routes to the 
critical points. 



The numerical results [11,12| agree well with an approximate analytic renormalization 



results obtained by Mao and Greene |jT3[ , except for the zero-coupling case in which the two 



area-preserving maps become uncoupled. Using an approximate renormalization method 
including truncation, they found three relevant eigenvalues, ^i = 8.9474, 62 = —4.4510 and 
63 = 1.8762 for the zero-coupling case |]14[. However they believed that the third one ^3 



is an artifact of the truncation, because only two relevant eigenvalues 61 and 62 could be 
indentified with the scaling factors numerically found. 

In this Brief Report we numerically study the critical behavior at the zero-coupling point 
in two symmetrically coupled area-preserving maps and resolve the discrepancy between the 
numerical results on the scaling of the coupling parameter and the approximate renormaliza- 



tion results for the zero-coupling case. In order to see the fine structure of period doublings, 
we extend the simple one-term scaling law to a two-term scaling law. Thus we find a new 
scaling factor ^3 = 1.8505 . . . associated with coupling, in addition to the previously known 
coupling scaling factor 82 = —4.4038 .... The numerical values of 62 and ^3 are close to 
the renormalization results of the relevant coupling eigenvalues 62 and 63. Consequently 
the fixed map governing the critical behavior at the zero-coupling point has two relevant 
coupling eigenvalues 62 and 63 associated with coupling perturbations, unlike the cases of 
other critical points. 

Consider a 4D volume-preserving map T consisting of two symmetrically coupled area- 



preserving Henon maps 11,12 



T: <^ 



(1) 



xi{t + 1) = -yi(t) + f{xi{t)) + g{xi{t), X2{t)), 

yi(t+l)=xi(t), 

X2(t+ 1) = -y2{t) + fix2{t)) + gix2{t),xi{t)), 

y2{t+l)=X2{t), 

where t denotes a discrete time, / is the nonlinear function of the uncoupled Henon's 
quadratic map [O, i.e.. 



/(^) 



1 — ax 



(2) 



and g{xi,X2) is a coupling function obeying a condition 



g{x,x) = for any x. 



(3) 



The two-coupled map (|T]) is called a symmetric map ||TT],|T2| because it is invariant under 
an exchange of coordinates such that xi ^^ X2 and yi ^^ y2- The set of all points, which are 
invariant under the exchange of coordinates, forms a symmetry plane on which Xi = X2 and 
yi = y2- An orbit is called an in-phase orbit if it lies on the symmetry plane, i.e., it satisfies 



xi{t) = X2(t) = x(t), yi(t) = y2{t) = y{t) for all t. 



(4) 



Otherwise it is called an out-of-phase orbit. Here we study only in-phase orbits. They can 
be easily found from the uncoupled Henon map because the coupling function g satisfies the 
condition (Rf). 



Stability analysis of an in-phase orbit can be conveniently carried out |TT|,|T2| in a set of 
new coordinates (Xi, Fi,X2, F2) defined by 



X2 



[xi +X2) 

2 
(xi - X2) 



Yi 



Yo 



2 

(2/1 - y2) 



(5a) 
(5b) 



2 ' 2 

Note that the in-phase orbit of the map (|l|) becomes the orbit of the new map (expressed 
in terms of new coordinates) with X2 = 1^2 = 0. Moreover the new coordinates Xi and Yi 
of the in-phase orbit also satisfy the the uncoupled Henon map. 

Linearizing the new map at an in-phase orbit point, we obtain the Jacobian matrix J 
which decomposes into two 2x2 matrices |Tl],|l2[: 



J 



Jo 



(6) 



Here is the 2x2 null matrix, and 



Ji 



J2 



( f (xo -1 ^ 



V 







/ 



( fix,) - 2G(Xi) -1 ^ 



V 



(7) 



(8) 



/ 



where f{X) = -^ and G(X) 



8g(Xi,X2) 
dX2 



Xi=X2=X 



. Hereafter the function G'(X) will be 



called the "reduced" coupling function of g{Xi, X2). Note also that the determinant of each 
2x2 matrix Ji {i = 1,2) is one, i.e., Det{Ji) = 1. Hence they are area-preserving maps. 

Stability of an in-phase orbit with period q is then determined from the g-product Mj of 
the 2x2 matrix J,-: 



q-l 

M, = l[J,{X,{t)), I 
t=o 



1,2. 



(9) 



Since Det{Mi) = 1, each matrix M, has a reciprocal pair of eigenvalues, Aj and Aj ^. Associate 
with a pair of eigenvalues (Aj, Aj~^) a stability index |jl6[, 

A = A, + A,-\ ^ = 1,2, (10) 

which is just the trace of Mj, i.e., pi = Tr{Mi). Since Mi is a real matrix, pj is always real. 
Note that the first stability index pi is just that for the case of the uncoupled Henon map 
and hence coupling affects only the second stability index p2- 

An in-phase orbit is stable only when the moduli of its stability indices are less than 
or equal to two, i.e., \pi\ < 2 for i = 1 and 2. A period-doubling (tangent) bifurcation 
occurs when each stability index pi decreases (increases) through —2 (2). Hence the stable 
region of the in-phase orbit in the parameter plane is bounded by four bifurcation lines 
associated with tangent and period-doubling bifurcations (i.e., those curves determined by 
the equations pi = ±2 for i = 0,1). When the stability index pi decreases through —2, the 
in-phase orbit loses its stability via in-phase period-doubling bifurcation and gives rise to 
the birth of the period-doubled in-phase orbit. Here we are interested in scaling behaviors 
of such in-phase period-doubling bifurcations. 

As an example we consider a linearly-coupled case in which the coupling function is 

g{xi,X2) = -{x2-Xi). (11) 



Here c is a coupling parameter. As previously observed in Refs. [Tl|JT^ , each "mother" 
stability region bifurcates into two "daughter" stability regions successively in the parameter 
plane. Thus the stable regions of in-phase orbits of period 2" {n = 0, 1, 2, ■ ■ ■) form a 
"bifurcation" tree in the parameter plane |]1^ . 



An infinite sequence of connected stablity branches (with increasing period) in the bifur- 
cation tree is called a bifurcation "route" |Tl|JT^ . Each bifurcation route can be represented 
by its address, which is an infinite sequence of two symbols (e.g., L and R). A "self-similar" 
bifurcation "path" in a bifurcation route is formed by following a sequence of parameters 
{an,Cn), at which the in-phase orbit of level n (period 2") has some given stability indices 



(pi,P2) (e.g., pi = —2 and p2 = 2) |TT|JT^ . All bifurcation paths within a bifurcation route 
converge to an accumulation point (a*, c*), where the value of a* is always the same as that 
of the accumulation point for the area-preserving case (i.e., a* = 4.136 166 803 904. . .), but 
the value of c* varies depending on the bifurcation routes. Thus each bifurcation route ends 
at a critical point {a*,c*) in the parameter plane. 

It has been numerically found that scaling behaviors near a critical point are characterized 
by two scaling factors, 6i and 62 [lllJl^. The value of ^i associated with scaling of the 



nonlinearity parameter is always the same as that of the scaling factor 6 (= 8.721 . . .) for 
the area-preserving case. However the values of 62 associated with scaling of the coupling 
parameter vary depending on the type of bifurcation routes. These numerical results agree 



well with analytic renormalization results [13|, except for the case of one specific bifurcation 



route, called the E route. The address of the E route is [(L, R, )°°] (= [L, R, L,R, . . .]) and 
it ends at the zero-coupling critical point (a*, 0). 

Using an approximate renormalization method including truncation, Mao and Greene 
P^ obtained three relevant eigenvalues, ^i = 8.9474, 62 = —4.4510, and ^3 = 1.8762 for 
the zero-coupling case; hereafter the two eigenvalues 62 and ^3 associated with coupling 
will be called the coupling eigenvalues (CE's). The two eigenvalues ^i and 62 are close to 
the numerical results of the nonlinearity-parameter scaling factor 6i{= 8.721 . . .) and the 
coupling-parameter scaling factor 62{= —4.403. . .) for the E route. However they believed 
that the second relevant CE 63 is an artifact of the truncation, because it could not be 
identified with anything obtained by a direct numerical method. 

In order to resolve the discrepancy between the numerical results and the renormalization 
results for the zero-coupling case, we numerically reexamine the scaling behavior associated 
with coupling. Extending the simple one-term scaling law to a two-term scaling law, we find 
a new scaling factor ^3 = 1.8505 . . . associated with coupling in addition to the previously 
found coupling scaling factor 62 = —4.4038. . ., as will be seen below. The values of these 
two coupling scaling factors are close to the renormalization results of the relevant CE's 62 
and ^3. 



We follow the in-phase orbits of period 2" up to level n = 14 in the E route and 
obtain a self-similar sequence of parameters {an,Cn), at which the pair of stability indices, 
(po,n,Pi,n), of the orbit of level n is (—2,2). The scalar sequences {a„} and {cn} converge 
geometrically to their limit values, a* and 0, respectively. In order to see their convergence, 
define 5„ = -—^ — and /i„ = —-^ — , where Aa„ = an — a„_i and Ac„ = c^ — c„_i. Then 

Attn Ac„ 

they converge to their limit values 6 and /i as n — > oo, respectively. Hence the two sequences 
{Aa„} and {Acn} obey one-term scaling laws asymptotically: 

Aa„ = C^'^^r", Ac„ = C('=V"" for large n, (12) 

where C^""^ and C^^^ are some constants, 5 = 8.721 ■ ■ -, and /i = —4.403 ■ ■ •. The values of 5 
and fi are close to the renormalization results of the first and second relevant eigenvalues 6i 
and 62, respectively. 

In order to take into account the effect of the second relevant CE S3 on the scaling of the 
sequence {Ac„}, we extend the simple one-term scaling law (p!2| ) to a two-term scaling law: 

Acn = Ci/Xj^" + 6*2/^2^" for large n, (13) 

where |/ii| > |/X2|- This is a kind of multiple scaling law |T^. Eq. (|r^) gives 



ACn = tiACn+1 - t2ACn+2, (14) 

where ti = /xi + fj,2 and ^2 = PiP2- Then /xi and fj,2 are solutions of the following quadratic 
equation, 

/i^ - hn + t2 = 0. (15) 

To evaluate fii and /i2, we first obtain ti and t2 from Ac„'s using Eq. (|I^: 

ACnACn+l — ACn-lACn+2 



Acl - ACn+lACn^i 

^c„_|_i — Ac„Ac„+2 



(16a) 
(16b) 



Note that Eqs. (|T3D-(|T6D hold only for large n. In fact the values of tj's and /ij's {i = 1,2) 
depend on the level n. Therefore we explicitly denote tj's and /Xj's by tj,n's and fJ.i^nS, 
respectively. Then each of them converges to a constant as n — > oo: 

lim tin = ti, lim /Xj „ = /Xj, i = 1, 2. (17) 

Three sequences {/ii,n}, {Ai2,n}, and {/i^ „//i2,n} are shown in Table |. The second column 
shows rapid convergence of /ii,„ to its limit values /ii (= —4.403 897805), which is close to 
the renormalization result of the first relevant CE (i.e., 62 = —4.4510). From the third and 
fourth columns, we also find that the second scaling factor fi2 is given by a product of two 
relevant CE's 62 and ^3, 

/i2 = |, (18) 

where 82 = yUi and 63 = 1.850 65 . It has been known that every scaling factor in the multiple- 
scaling expansion of a parameter is expressed by a product of the eigenvalues of a linearized 



renormalization operator |T3 . Note that the value of ^3 is close to the renormalization result 
of the second relevant CE (i.e., ^3 = 1.8762). 

We now study the coupling effect on the second stability index p2,n of the in-phase orbit of 
period 2" near the zero-coupling critical point (a*, 0). Figure |l| shows three plots of p2,n('3'*, c) 
versus c for ra = 4, 5, and 6. For c = 0, p2,n converges to a constant pg (= —2.543 510 20 . . .), 



called the critical stability index |jT2[, as n ^ 00. However, when c is non-zero p2,„ diverges 

dp2,n 



as n -^ 00, i.e., its slope Sn 



dc 

n -^ 00. 



) at the zero-coupling critical point diverges as 

(a*,0) 



The sequence {Sn} obeys a two-term scaling law, 

Sn = Div'^ + D2V2 for large n, (19) 

where |z/i| > \v2\- This equation gives 

Sn+2 = nSn+l -r2Sn, (20) 



where ri = ui + z/2 and r2 = z^it'2- As in the scahng for the couphng parameter, we first 
obtain ri and r2 of level n from SnS: 

C C C C q2 Q Q 

'^l,n — ^ 7; 7; ; ^2,n — '^ 7; ^; • {■^^J 

Then the scaling factors ui^n and z/2,n of level n are given by the roots of the quadratic 
equation, z/^ — ri ,ii/„ + r2,n = 0. They are listed in Table |T| and converge to constants z/i 
(= —4.403 897805 09) and z/2 (= 1.850 535) as n -^ 00, whose accuracies are higher than 
those of the coupling-parameter scaling factors. Note that the values of z/i and z/2 are also 
close to the renormalization results of the two relevant CE's 62 and ^3. 

We have also studied several other coupling cases with the coupling function, g(xi,X2) = 
-{x2 — x") (n is a positive integer). In all cases studied (n = 2,3,4,5), the scaling factors 
of both the coupling parameter c and the slope of the second stability index p2 are found to 
be the same as those for the above linearly-coupled case (n = 1) within numerical accuracy. 
Hence universality also seems to be well obeyed. 
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TABLES 
TABLE L Scaling factors /ii^„ and ^2,n in the two-term scaling for the coupling parameter are 

Ml n 

shown in the second and third columns, respectively. A product of them, — '—, is shown in the 



fourth column. 








n 


/^l,n 


M2,n 


/"2,n 



5 -4.403 908128 10.437 4 1.85817 

6 -4.403 899 694 10.465 9 1.853 09 

7 -4.403 898 736 10.458 2 1.854 46 

8 -4.403 897 867 10.474 8 1.85152 

9 -4.403 897 847 10.473 9 1.85168 

10 -4.403 897 806 10.478 4 1.850 89 

11 -4.403 897 807 10.478 6 1.850 85 

12 -4.403 897 805 10.479 7 1.850 65 

TABLE IL Scaling factors ui^n and V2,n in the two-term scaling for the slope of the second 
stability index are shown. 

n i^l,n l^2,n 

5 -4.403 898453 59 1.851433 5 

6 -4.403 897 730 29 1.850 782 6 

7 -4.403 897813 85 1.850 603 6 

8 -4.403 897804 07 1.850 553 8 

9 -4.403 897805 21 1.850 540 

10 -4.403 897805 07 1.850 5361 

11 -4.403 897805 09 1.850 535 

12 -4.403 897805 09 1.850 534 9 
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FIGURES 
FIG. 1. Plots of the second stability index p2,n{o* ^ c) versus c for n = 4, 5, 6. 
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